Mon. Not. R. Astron. Soc. OOP. [TlfTSl C2008') Printed 16 March 2009 (MN WT^ style file v2.2) 



On the peculiar momentum of baryons after Reionization 



Carlos Hernandez-Monteagudo^* and Shirley Ho^ 

^ Max Planck Institut fiir Astrophysik (MPA), Karl Schwarzschild Str.l, Garching bei Miinchen, D-85741, Germany 
^Lawrence Berkeley National Laboratory, Berkeley, CA 94704, USA 



ABSTRACT 

The peculiar motion of ionized baryons is known to introduce temperature 
anisotropies in the Cosmic Microwave Background radiation (CMB) by means of the 
kinetic Sunyaev-Zel'dovich effect (kSZ). In this work, we present an all sky compu- 
tation of angular power spectrum of the temperature anisotropies introduced by kSZ 
momentum of all baryons in the Universe during and after reionization. In an at- 
tempt to study the bulk flows of the missing baryons not yet detected, we address 
separately the contribution from all baryons in the Inter Galactic Medium (IGM) and 
those baryons located in collapsed structures like groups and clusters of galaxies. In 
the first case, our approach provides a complete, all sky computation of the kSZ in 
second order of cosmological perturbation theory (also known as the Ostriker-Vishniac 
effect, OV). Most of the power of OV is generated during reionization, although it has 
a non-negligible tail at low redshifts, when the bulk of the kSZ peculiar momentum 
of the halo (cluster -I- group) population arises. If gas outside halos is comoving with 
clusters as the theory predicts, then the signature of the bulk flows of the missing 
baryons should be recovered by a cross-correlation analysis of future CMB data sets 
with kSZ estimates in clusters of galaxies. For an ACT or SPT type of CMB exper- 
iment, all sky kSZ estimates of all clusters above 2 x 10^^ h~^MQ should provide a 
detection of dark flows with signal to noise ratio (S/N) of ^ 10, (S/N ~ 2.5 — 5 for 
2,000 - 10,000 square degrees). Improving kSZ estimates with data from Large Scale 
Structure surveys should enable a deeper confrontation of the theoretical predictions 
for bulk flows with observations. The combination of future CMB and optical data 
should shed light on the dark flows of the nearby, so far undetected, diffuse baryons. 
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1 INTRODUCTION 

The standard cosmological model predicts the statistical 
properties and evolution of the large scale matter distribu- 
tion. Initially small inhomogeneities in the matter density 
field grow under gravity in an expanding universe, giving rise 
to a cosmic web made of filaments, sheets and superclusters 
of galaxies that we see today. Large, comoving flows of mat- 
ter fall onto deep potential wells under gravity, with typical 
speeds of ~ 150 km s~^ in th e standard ACDM cosmological 
model (|Hinshaw et al.ir2008l ). These peculiar velocities ( pe- 
culiar with respect to the Hubble flow induced by the univer- 
sal expansion) cannot be easily measured. I ndeed, despite re- 
cent c laims pending for co nfirmation fe.g:.. lKashlinskv et al.l 
l2008l : IWatkins et al.ll2008l '). there is no conclusive observa- 
tional evidence for such bulk fiows in the present universe. 



Nevertheless, the Cosmic Microwave Background 
(CMB) radiation has provided an accurate measurement 
of these bulk fiows when the universe was 380,000 years 
old (2 ~ 1,050). Before most of the electrons recombined 
with the protons and the CMB could propagate freely, the 
electron plasma and the CMB radiation were in thermal 
equilibrium since Thomson scattering kept both fluids cou- 
pled. The free electrons Thomson scattered CMB photons, 
leaving Doppler imprint according to their relative veloci- 
ties caused by the oscillation of the plasma-photon fluid in 
the gravitational poten tial wells l|Sunvaev fc Zeldovicbll970l : 
iHu fc Sugivaiiial Il995l ). Since energy and density fluctua- 
tions were very small at that epoch, the linear theory is able 
to describe very accurately the anisotropy pattern that was 
introduced by Thompson scattering on the CMB, and this 



has been confirmed by observations (e.g. Mausko pf et al 
2OOOI : iMiUer et al1l2002l : ICrainge et "alll2003l : iHinshaw et al 
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At redshift 2: ~ 10 the first stars began to ionize hydro- 
gen in the Intergalactic Medium (IGM), and therefore CMB 
photons interacted again with the free electrons via Thom- 
son scattering. In this scenario, electron peculiar velocities 
introduced new temperature anisotropies on the CMB. Like- 
wise, due to the anisotropic nature of Thomson scattering, 
the CMB was also linearly polarized on the large angular 
scales, as it has also b een confirmed by the WMAP experi- 
ment, ijHinshaw et al ] [2008. ). These epochs describe the last 
scenario where all baryons have been observed: at more re- 
cent epochs, the amount of baryon matter that we are able 
to detect is at most of half of w hat is observed during bot h 
recombination and reionization (|Fukugita fc Peebles! l20o3 ). 
This constitutes the missing baryon problem. 

In this work, and jus t as in 

iHernandez-Monteaeudo fc Sunvae we investi- 

gate bulk flows and the missing baryons utilizing Thomson 
scatterings between CMB photons and free electrons at 
low redshift. This mechanism is refe rred to as the kinetic 
Suny aev-Zel'dovich effect (kSZ, ISunvaev fc ZeldovichI 
1 19721 ). ii it occurs in galaxy cl usters, but it is called the 
Ostriker Vishniac effect (OV, lOstriker fc Vishniad Il986l l 
in the context of cosmological second order perturba- 
tion theory, (for which the electron density and velocity 
are continuous fields). We shall follow this terminology, 
although in some occasions we may use either when 
referring to both effects for simplicity. We shall restrict 
ourselves to the impact of Thomson scattering on the 
intensity or temperature anisotropies of the CMB, since 
the effect on the CMB polarization i s alr eady studied 
in iHernandez-Monteagudo fc SunvaevI (l2P08) . Unhke in 
iHernandez-Monteagudo et al.l l|2006l '). we shall not consider 
the spectral distortions that the hot fraction of the missing 
baryons introduce in the CMB black body spectrum. 
Previous eff orts have been made i n the computation of the 
OV effec t (lOstriker fc Vishniad 1 19861 : iDodelson fc JubasI 
I1995I : IHuI 120001 : ICastrdl2003h . some of them in the context 
of the non-Gaussian signal they give rise to. In some 
recent works, the kSZ in halos has been used as a probe 
for Dark Energv ( | Bhattacharva fc Kosowskvi i200S ) , while 
iMoodlev et al.l 1 20081 ) hav e used the gas in groups to probe 
the missing baryons, and iLed (|2009l ) studies the prospects 
of constraining patchy reionization. Here, we present the 
first full sky projection of the temperature anisotropies 
induced by peculiar momenta of ionized baryons during 
and after reionization up to the present: if applied to the 
smooth electron distribution our calculations provide a full 
description of the OV effect, while when applied to the 
halo distribution they describe the contribution of galaxy 
group/cluster correlation on the kSZ effect. The formalism 
also permits the computation of the cross correlation 
between the OV and the kSZ effects. At the same time, our 
computations provide the linear approximation under which 
the electron density is approximated by its average value. 
We can then compare our code with s tandard linear theory 
Boltz mann code s such as CMBFA ST (jSeliak fc Zaldarriagal 
1 1996! ) or CAMB (jLewis et al.!!2000 h. as well as normalize our 
predictions to the linear fluctuation amplitudes measured 
by WMAP. Finally, we analyze our results in the context 
of future CMB and Large Scale Structure observations, 
and provide predictions for the detectability of bulk flows 
and missing baryons in the light of those future data 



sets. The paper is structured as follows: in Section [2] we 
describe the line of sight projection of the baryon peculiar 
momentum, and this is discussed in Section|3]in the context 
of its application to the halo population. The amplitude 
of the resulting power spectra are given in Section |4j and 
their interpretation in terms of bulk flow/missing baryon 
detection is provided in Section [S] We discuss our main 
results in Section [6] and conclude in Section [7] 



2 THE SECOND ORDER KSZ EFFECT (OV 
EFFECT) 

The computation of the OV/kSZ effect requires the projec- 
tion along the line of sight the product of electron peculiar 
velocity and density. If we consider the spatial variations of 
the electron density, then this product becomes a second or- 
der quantity in cosmological perturbation theory provided 
that peculiar velocities are very small (u/c <^ 1). I.e., the 
CMB temperature anisotropies introduced by the peculiar 
motion of the scatterers reads 



5T 



with tt the Thomson optical depth 



drj a{rj)aTne{ri), 



(1) 



(2) 



a{ri) = 1/(1 -I- 2(77)) the cosmological scale factor, or 
the Thomson cross section and n^{ri) = ne(l -I- 5^) is 
the electron number density in terms of its average value 
ijie) and its spatial fluctuations (5e). The symbol r] de- 
no tes conformal time or, equiva l ently, comoving distance. 
In !Hernandez- Monteagudo et al" I (120061 ) (hereafter HVJS) a 
full sky description of the kSZ generated by the low redshift 
galaxy cluster populat i on wa s provided. However, as noted 
later bv ! Gordon et al.! (|2007l ). those computations observed 
only the diagonal part of the Fourier peculiar velocity tensor 
(iij^iiq*) (where k and q are Fourier wave vectors and i and 
j refer to one of the three spatial components of the peculiar 
velocity vectors), neglecting all non-diagonal components. In 
Appendix A, we present a full sky computation considering 
all terms in this tensor, and find that the complete treat- 
ment introduces small changes to the results given in HVJS. 
The results of the Appendix are summarized here. 

In brief, the second order statistic for the electron 
momentum (pe oc UeV^) is proportional to the quantity 
{ne,iVe^ine,2Ve,2) ■ If We,i is coustaut, thcu it simplifies to 
n'i{ve,iVe,2) , and this is the vv term. Furthermore, it is pos- 
sible that Ue is driven by Poissonian statistics, in which case 
{nc,iVe^ine^2Ve,2) OC ne{?;e,iWe,2) (Poissou term). If however 
Tie,! is correlated to the velocity field, then one must use the 
cumulant expansion theorem, and three new terms arise, of 
which one is positive, another one is zero and the third one 
is negative. In all cases, the correlation function C(ni,n2) 
for each term can be written as Fourier integrals over some 
wave vector k of functions with an angular dependence of 
either the type (k • ni)(k • £12) or (fii • 112). Those contribu- 
tions with (k-ni)(k-n2) type dependence will give negligible 
contribution at the small scales, i.e., the corresponding an- 
gular power spectrum multipoles Ci will go to zero as the 
multipole I grows. In this limit, only the contributions with 
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the (fii • 112) type angular dependence will survive. Let us 
consider the above more carefully in case by case basis. 

2.1 The linear terms 

If the electron density is constant, then kSZ/OV effect is 
linear and the corresponding power spectrum will contribute 

to the total with a term that wc shall name as vv. The 
amplitude of the vv term is given by 
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(2Z + 1)2- 
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In this equation, the terms Tf^ correspond to the projection 
on multipolc I of the line of sight integrals of the rodshift do- 
pendent source term that includes the growth factor of the 
peculiar velocities "D^ (z) and the average Thomson visibility 
function tt exp (— rr) (with t't ~ CTfic a{z) the Thomson 
opacity, dependent on the Thomson cross section, the av- 
erage electron density (rie) and the scale factor a{z)). The 
label i refers to different lines of sight nj, i = 1, 2. The linear 
peculiar velocity growth factor Vy can be expressed in terms 
of the linear density growth factor (I'd), 



X)„ = H{z] 



dVs 



dz 



and T>s[z) is expressed in terms of the integral 
dz 



Vs{z) oc H{z) 



(4) 



(5) 



with H{z) the Hubble function and Ds normalized to Ds{z = 
0) = 1 .The spatial clustering properties of the peculiar ve- 
locity field are encoded in the integral of the velocity power 
spectrum, which is proportional to P{k)/k'^ {P{k) is linear 
matter power spectrum) . Provided that all lines of sight are 
statistically equivalent and for high I there exists an asymp- 
totic value for the spherical Bessel functions ji(x), it is easy 
to find that at the small scales the vv term vanishes, i.e., 

lim Cr = 0. (6) 

This term is linear, and is therefore present in standard 
Boltzmann codes like CMBFAST or CAME. 

In the case of halo induced kSZ, the discrete nature of 
halos introduces a different type of anisotropy. The actual 
amplitude of the kSZ is proportional to the number of ha- 
los present along the line of sight, therefore kSZ fluctuations 
obey Poisson statistics in this case. The variance of the num- 
ber of objects is proportional to its mean, and so the angular 
power spectrum is also linear in the average halo number 
density. The Poisson term has two contributions, Cf^ and 



Ci^A = (^) j k^dkji{k) X 



2Z-hl 



-'•i-l,l-'-i-l,2 



(7) 



and it does not vanish on the small scales when both lines 
of sight are statistically equivalent, although it will drop to 
zero in the angular range that goes well beyond the typical 
halo profile (halos are assumed to have no substructure. 

The function Ji(k) is an angular integral involving the 
window functions from the Fourier transform of the halo pro- 
file and the matter power spectrum. The terms Xfi express 
the multipole integral along the line of sight the visibihty 
function and the velocity growth. They are proportional to 
\/rih {rih is the halo rmmbcr density), so that their product 
is linear in the halo density, as Poisson statistics requires. 
On the other hand, the term Cfg reads 



Cf,B 
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k'dk J2{k)x 



T'- T'- (Z -I- \)l p 
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jl + 1)^ T-p , I' 

~ (2[TiF ■ '(21 + 

This term vanishes again at high l-s, 
lim = 0, 



2^-1,22^-1,2 1 ■ (8) 



(9) 



so in this multipole regime Cf CC^A- The function J2{k) 
is another A; dependent function similar to Ji{k). 



2.2 The second order terms 

Both the vv and the Poisson terms above are linear, and 
hence their angular power spectra are proportional to the 
linear scalar power spectrum (although wc prefer to write 
them in terms of the linear matter power spectrum). Sec- 
ond order terms involve the product of density and peculiar 
velocity, and thus resulting angular power spectra contain 
terms that can be interpreted as convolutions of the initial 
scalar power spectrum. In the Appendix, we demonstrate 
that three different terms of this nature arise, namely the 
vidi — vidi term, the t;it;2 — didi term and the vidi — vidi 
term. (The integer subscript here denotes one of the two 
lines of sight). The first one turns out to give zero contribu- 
tion, so for simplicity we shall refer to the latter two as the 
vv — dd and vd — vd terms, respectively. 

The vv — dd term comes from two contributions. 



cr 
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and, as the vv term above, it vanishes on the small scales. 
The expression for C^-'^'' is 

c;;;b-'"'= (^) / k^dkB{k) x 
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21 + 1 



(12) 



and dominates the contribution to (Ji^^'^'^ at high multi- 
poles. Contrary to Ji{k), J'2{k) above, the functions A{k) 
and B{k) are angular integrals involving the convolution of 
two matter power spectra, apart from the window functions 
from Fourier transform of the halo profile. The terms X"""'^'^ 
provide the line of sight integral on multipole / of the visi- 
bility function, the matter bias factor (matter bias factor is 
set to unity unless when we are considering halos) and the 
growth factors of density and velocities, respectively. 



The term C, 



is very similar to C" 



upon two different contributions, 

Cvd — vd y^vd — vd . y^vd—vd 
I — '^l A + L"; R 1 



it is built 



(13) 



and the first one also vanishes on the small scales, as can be 
see from this expression when we go to large Z: 
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The second contribution is given by 
2 
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21 + 1 



(15) 



The functions C{k) and 2?(fc) are again angular integrals 
which contain the convolution of two power spectra, spher- 
ical Bessel functions and window functions associated with 
the size of the regions where peculiar velocities are mea- 
sured. 



3 GAS IN HALOS VERSUS GAS IN THE IGM 

Before we show the amplitudes of the OV/kSZ power spectra 
presented in the previous section, we compare the effect of 
Thomson scattering as it arises in collapsed halos to the OV 
generated in the intergalactic medium (IGM). It is relevant 
to understand under which conditions would collapsed gas in 
halos contribute more in terms of CMB angular anisotropics 
than the gas in the IGM. In both scenarios, Thomson scat- 
tering introduces anisotropics according to 



ST 
To 



(n) = dr]TT 



(16) 



where we have approximated the Thomson visibility func- 
tion A = rTexp(— tt) by the opacity, A ~ fx, since tt 
is ~ 0.09 in the redshift regime of interest (jHinshaw et al.l 
l2008l ). 

We assume that the peculiar velocity of gas with respect 
to the CMB is very similar for both gas in IGM and halos, 
the comparison of the OV and kSZ effects must then be 
focused on the optical depths, i.e., in the line of sight integral 




0.10 1.00 
Redshift z 

Figure 1. Top: Opacities in a WMAP V Universe due to a 
smooth distribution of electrons (solid line), all halos more mas- 
sive than 5 X 10^h~^MQ (dotted line) and all halos more mas- 
sive than 5 X 10^^ Mq (dashed line). Bottom: Thin lines re- 
fer to the M > 5 X 10^/i~^Mq halo population, thick lines to 
M > 5 X 10^^ Mq. Dotted lines: Bias factor evolution versus 
redshift. Solid lines: Volume fraction filled by each halo popu- 
lation. Dashed lines: Product of volume fraction times electron 
overdensity (with respect to background). 



of t't in each case. As shown in HVJS, the effective opacity 
generated in halos at a given epoch is given by 



{TT)halo{z) = 



dM ■^TT,halo{M,z)V(A4,z) = 



dM iT,halo{M,z)f^ol(M,z) 



(17) 



where dn/dM is the halo mass function (w e choose 
Sheth and Tormen (ST. ISheth fc TormenI [l999h . V[M,z) 
is the proper volume occupied by gas in those halos, and 
TT,haio{M, z) is the opacity generated by a halo of mass M 
at redshift z. The effect is in the end sensitive to the to- 
tal number of electrons in a given halo population, and can 
be expressed by the average halo optical depth times the 
volume fraction of halos. 

The halo population opacity can be compared directly 
to the opacity associated to the smooth distribution of all 
electrons. The top panel of Figure ([1} displays the opacity in 
three scenarios: (i) a, smooth distribution of electrons (solid 
line)Q, (ii) electrons within all halos of mass greater than 
5 X 10^^ Mq (dashed line), and (Hi) electrons within all 
halos of mass greater than 5 x 10^ Mq (dotted line). The 
most massive halo population provides a very low opacity 
due to the small amount of cosmological volume it takes. 
The latter halo population contains almost all the electrons 
in the Universe, and therefore its opacity is very close to the 
one in case (ijj. 

^ We assume a sudden reionization model within WMAP V cos- 
mology, for which tt = 0.084 ± 0.016 and Zreio = 10.8 ± 1.4 
2 According to the halo model, all mass in the Universe must be 
found in collapsed halos 
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We have ignored the effect of the bias until now. Dark 
matter halos can either be biased or anti-biased tracers of 
the underlying matter distribution, and we need to include 
the bias factor when we compute the angular anisotropies. 
In the bottom panel of Figure ([l| the bias factors for the 
two halo populations considered (M > 5 x 10^ Mq and 
M > 5 X 10^^ /i~^M0) are shown by the dotted thin and 
dotted thick hnes, respectively. Th ese bias factors were com- 
puted by using the expressions of IMo fc Whit3 (|l996l ). and 
integrating them with the ST mass function: 

K^)^J'^M^KM,.)/ jdM^. (18) 

Although at high redshifts bias factors grow above unity, the 
volume fraction occupied by the halo population are very 
low: this is shown by bottom solid lines, (thin line for M > 
5 X lO^/i"^M0, thick line for M > 5 x 10"/i"^A/o). The 
volume fraction is so low that it cannot be compensated by 
gas overdensity within halos: the product of gas overdensity 
and the volume fraction is below unity at all redshifts (as 
shown by the thin (M > 5 x IO^/i'^Mq) and thick (M > 
5 X lO^^/i-^M©) dashed lines). The thin solid green line 
shows the level of unity. 

According to these results, it seems that the anisotropy 
due to halos will always be below the level of anisotropy 
induced by a smooth and continuous distribution of elec- 
trons. There are however two caveats: (i) Halos might show 
bias in their peculiar velocities, although this i s unlikely t o 
change the picture by more than a 30%-40% |Pee? 2006), 
(ii) on the small scales, halos follow Poissonian statistics, 
and this will add a considerable amount of anisotropy, as 
we s hall see below. Th ese results differ from previous works 
(e.g., IHuI|2000I : ICastrd l2003) where the power spectrum for 
kSZ in non linear structures had been computed. In previous 
works, the matter power spectrum was simply replaced by 
the non linear power spectrum, but no changes were intro- 
duced in the related opacity. No attention was paid to the 
Poisson term either, which, in the case of halos, turns out 
to be dominant. This difference must be motivated by the 
fact that, they were computing the combined effect of both 
linear and non-linear density perturbations at high redshift. 
Here we try to separate the OV contribution from the one 
generated in halos, and the contribution to the latter arises 
mostly at low [z < 2) redshift. 



4 THE POWER SPECTRUM OF THE OV/KSZ 
EFFECT 

4.1 Normalization to WMAP observations 

The WMAP satellite has provided all sky measurements of 
the temperature and polarization anisotropies of the CMB 
ijHinshaw et al.ir2008l l. This fluctuation field is so small that 
should provide a normalization for the linear level of de- 
parture from perfect anisotropy (higher order contributions 
should be negligible) . On the large scales, the theory predicts 
that the CMB temperature anisotropies must caused mainly 
by two physical phenomena: gravitational redshift/blueshift 
of CMB photons due to the passage through different (and 
eventually evolving) gravitational potentials, and Thomson 
scattering induced by free electrons during and after reion- 
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Figure 2. Comparison of the WMAPV-normalized prediction of 
the vv term by CMBFAST (solid line) with our computation (red 
filled circles). 



ization. This theory predicts the relative weight of these two 
contributions, and WMAP measurements constitute accu- 
rate measurements of the combined effect of these two contri- 
butions. Therefore, if the theory is correct, by using WMAP 
observations it is possible to estimate the amplitude of each 
of the two contributions dominating at the large angular 
scales. In particular, WMAP data can be used to normalize 
the level of anisotropies generated by Thomson scattering 
during and after reionization. 

In Figure ([2|), the black solid line displays the vv term 
which accounts for the Thomson scattering contribution to 
the low / CMB anisotropy, (see Equation (|3}). This has been 
produced by feeding the cosmological parameter set given for 
WMAPV-hSNALL-hBAO0 to a modified version of CMB- 
FAST. This code provides both the total amplitude of the 
angular power spectrum and the C"" term. The total power 
spectrum was normalized to fit WMAPV observations, and 
the vv term was scaled accordingly. As shown by the solid 
line in Figure ([2}, the maximum of the vv angular band 
power spectrum is reached at I ~ 20 — 30, with an amplitude 
close to ~ 22 (fiK)^ . Note that here power spectrum ampli- 
tudes are given in terms of quantity = 1(1 + l)Ci/(2n). 
The filled red circles display the result of our computation 
of Equation ([Sj: the accurate computation of this term is 
particularly difficult, since at large Vs it drops to zero as 
the vanishing difference of several integrals. Small numer- 
ical errors in the computation of these integrals introduce 
some scatter with respect to the CMBFAST prediction, but 
only at moderate and high I's. For I < 100, the accuracy of 
our code is better than 8%, and we should expect a simi- 
lar or better level of accuracy for all other terms which do 
not vanish as the vv term, (we remark that all other terms 
become dominated at high Ps by at least a non- vanishing 
contribution). Note that this comparison allows us normal- 
ize the output of our code not only for the linear predic- 

^ See jhttp : //lambda . gsf c ■ nasa ■ gov/product/map/dr3/params | 
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tions,but also for the second order predictions of vv — dd 
and vd — vd. Therefore, with the exception of the percent 
level uncertainty of our integrals, all sources of uncertainties 
are associated to the modeling of the mass function, the bias 
factor and the halo profile description. 



velocities which would boost the kSZ amplitudes in halos. 
However, to this point, there is very little observational ev- 
idence for these non-linear component (as well as for the 
linear one). 



4.2 Results 

The amplitude of the OV/kSZ power spectra in different 
scenarios are displayed in Figure (|3]). Note that T)^^^ = 
yJl{l + l)Cl/{2-K) in ^iK. The left panel considers all elec- 
trons in a diffuse and continuous phase, and should be a fair 
description of the electron distribution during reionization. 
This case adopts the opacity displayed by the thick solid 
line in the top panel of Figure ([ij , for which the anisotropy 
is mostly generated at z ~ 6 — 10. The linear vv term we 
use to calibrate the amplitude of the fluctuations peaks at 
~ 4 — 5 for Z ~ 10 — 20, but then drops rapidly for 
I > 100. On the contrary, the sum of the vv — dd and vd — vd 
terms (which corresponds to what is commonly understood 
as the OV eff'ect) becomes dominant at I ~ 200, and shows a 
rather fiat pattern on the small scales. The OV term crosses 
the linear power spectrum generated during recombination 
at Z ~ 4, 000, with an amplitude between 1 — 2 ^iK. Provided 
that there is no way to spectrally distinguish the power due 
to OV during reionization from the power originated at the 
last scattering surface, the crossover of the OV over the lin- 
ear power spectrum is particularly significant, since it pro- 
vides a direct observable of the reionization epoch. Its am- 
plitude is within the nominal sensitivity of third generation 
CMB experim ents hke ACT llFowler fc ACT CollaborationI 
l2006l ) or SPT (|Ruhl et al.ll20M V 

Indeed, the amplitude of the OV effect is considerably 
higher than the amplitude corresponding to the vv — dd 
and vd — vd terms generated by halos. It is well above 
the anisotropy level generated by the galaxy cluster pop- 
ulation which experiments like ACT and SPT are targeting 
(M > 2 X 10" h-^MQ, middle panel in Figure (O). Even 
for M > 5 X 10^ h~^MQ, the second order kSZ terms lie 
about one order of magnitude below the OV contribution at 
Z ~ 3, 000 (which translates into two order of magnitudes in 
the power spectrum) . This proves that most of the power of 
the OV is actually generated at high redshifts, since the OV 
contribution from low redshifts is well sampled by the kSZ 
corresponding to masses A/ > 5 x 10^ h~^MQ (note that 
the opacities of the OV and the kSZ for this halo popula- 
tion are very similar at low z, say, at z ~< 1, see Figure 
([1])). The difi'erence in the amplitudes of the OV and halo 
kSZ contributions can therefore be assigned to z > 1. 

However, by looking at the right panel of Figure @ it 
becomes clear that in order to distinguish the high-z OV con- 
tribution from the low-z kSZ halo contribution one must take 
into account the Poisson term. This can be achieved by mask- 
ing out all massive, clearly detected halos: after removing the 
Poisson contribution from clusters above 2 x 10^* h~^MQ, 
the Poisson term shown in the right panel of Figure ^ drop 
by a ~ 40%, which, according to our computations, should 
be below the OV level. We must note that the halo pecu- 
liar velocities assumed here are close to the linear limit (we 
have introduced only a mild bias of 30%), but reality might 
be very different due to the presence of thermal, non linear 



5 PROSPECTS FOR DETECTING BULK 
FLOWS/MISSING BARYONS IN FUTURE 
CMB DATA 

Of all kSZ terms shown in Figure ((3)l, one of the easiest to 
observe is the Poisson term at i ~ 2, 000 — 3, 000. Its am- 
plitude is almost as high as the combination of the vv — dd 
and vd — vd terms from the OV and more importantly, it 
is associated with massive halos that can be identified at 
other frequencies. Nevertheless, one must note that measur- 
ing the kSZ Poisson term at positions of galaxy clusters does 
not properly test the theoretical predictions of the model on 
the bulk flows. This term is sensitive to the average velocity 
dispersion of those objects (which is likely to be contami- 
nated by non-linear contributions) but tells nothing about 
the spatial correlations of the peculiar velocity fleld. 

In the context of the standard model, it was shown in 
HVJS that typical correlation lengths for bulk flows are 
about 20 - 40 h'^Mpc in comoving units. This means 
that all electrons within a sphere of radius from 20-40 
/i~^Mpc share a similar peculiar velocity. If this is an ac- 
curate description of reality, then one can use the kSZ de- 
tected in the most massive halos as a template for the pe- 
culiar velocity field in their surroundings. This idea was 
used in iHernandez-M onteagudo fc SunvaevI (|2008l ) to track 
the missing baryons by kSZ - E polarization mode cross- 
correlation, and can also be applied in this context. Let us 
assume that future high resolution CMB experiments pro- 
vide accurate measurements of the kSZ in halos more mas- 
sive than, say, 2 x 10^** h~^MQ, so that a template of the 
kSZ on the sky (hereafter Mfesz) can be built. If this tem- 
plate is then cross-correlated with the CMB sky (hereafter 
TcA/s), then not only does signal come from the known 
clusters, but also from the surrounding gas. Fundamental 
sources of noise for this cross-correlation are, to first order, 
the primordial CMB fluctuations and the instrumental noise 
(which usually becomes relevant on the scales of the PSF). 
This cross-correlation can be picked-up via a matched filter 
(MF) approach, 



- CMB 



C-^Mfesz 



MiszC-^Mksz 



(19) 



where the CMB signal is assumed to intrinsically contain 
some kSZ component proportional to the kSZ template, i.e., 
a'M.kSz- The covariance matrix C contains the noise sources 
(namely primordial CMB and instrumental noise). The sig- 
nal to noise ratio (S/N) of this cross correlation is given 
by 



S_ 

N 



(20) 



If the kSZ template is written in multipole space, then the 
last equation on the S/N can be written for each multi- 
pole I. In this basis, the covariance matrix C is diagonal for 
white instrumental noise A^;, (C);,;/ = {Cf^''^ +Ni)5i^i, , and 
Equation ((20} reads 
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1st and 2nd [OV] orders kSZ in halos M > 2x10" M,/h kSZ In holds M > 5x10' M,/h 




Uultipole t Multlpole I Multlpole I 



Figure 3. Angular power spectra of the OV/kSZ effect in fiK units, {T>y^ = l{l + l)Ci / {2it)). The red filled circles provide the 
difference in amplitudes between the vv — dd and the vd — vd term, since the latter is negative. Left panel: Anisotropy amplitude generated 
by a smooth distribution of electrons during and after reionization (OV effect). Middle panel: First and second order kSZ amplitude 
generated a cluster population with M > 2 X 10^* h~^MQ. Right panel: First and second order kSZ amplitude generated by a cluster 
population with Af > 5 X 10^ H-'^Mq. 




Figure 4. S/N achieved by cross correlating a kSZ/peculiar velocity template with an all sky CMB map. The kSZ template is built 
with estimates of peculiar motions in all clusters more massive than 2 X 10^'^h~^ Mq. Dashed line refers to the case in which cluster 
are included in the cross-correlation analysis, whereas solid lines display the case where known halos are masked on the CMB maps, 
so that only gas outside clusters contributes to the cross-correlation. Black line refers to an ACT-like experiment, red line refers to 
an experiment with the sensitivity and angular resolution of Planck's 217 GHz channel. Filled triangles (squares) refer to an ACT-like 
experiment covering 10,000 (2,000) square degrees. The three panels refer to different survey depths: Zmax = 0.7 (left), Zmax = 0.9 
(center), and all z (right). 



For non-full sky coverage {fsky < 1), there are not (21+1) de- 
grees of freedom per multlpole I, but roughly (21 + 1) fsky In 

1/2 

this case, the last equation should be multiplied by f^j.^ . The 
term C'j^'^^z accounts for the halo-halo kSZ auto-correlation 



(i.e., the Poisson, vv, vv — dd and vd — vd terms), whereas 
Cffcsz expresses the cross-correlation between the kSZ in- 
duced in halos and the OV generated by the continuous 
electron distribution on the large scales (i.e., mostly out- 
side halos). This cross-power spectrum can be easily com- 
puted by introducing, in our angular cross-correlation func- 
tion C(ni,n2), the halo visibility function in the los inte- 
grals Ii-s along direction fii, and the OV visibility function 
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along direction n2, (see Appendix). It is worth to stress that 
in this case there is no contribution from the Poisson term, 
since this term appears exclusively for the halo auto cor- 
relation. If we are searching for the kSZ signal of the gas 
outside halos, one should mask the known clusters, which is 
equivalent to dropping the C^'^gz term in the numerator of 
Equation ((21}. 

The result of adding all the S /N below a given multipole 
I is shown in Figure (jl]). Black lines refer to an ACT-SPT like 
all-sky CMB experiment with a sensitivity level of ~ l^K 
per resolution element of ~ 1 arcmin. The red color lines 
(both solid and dashed) display results for (all sky, fsky ~ 
1) Planck's HFI 217 GHz channel. If kSZ clusters are not 
masked, then the cross-correlation test yields a very high 
S/N, as displayed by the dashed lines. Instead, if known kSZ 
clusters are masked, then the S/N drops to a level than only 
an ACT-SPT type CMB experiment would be able to detect 
it, since it mostly comes from the small angular scales. This 
is shown by the solid lines, which depict the ideal case of 
full sky coverage {fsky = 1). On the contrary, filled triangles 
and squares refer to an ACT-SPT type experiment covering 
2,000 and 10,000 square degrees, respectively. We see that 
these scenarios lie in the limit of detectability. 

The latter case describes the situation in which kSZ 
measurements in galaxy clusters are used to build the kSZ 
template and search for the baryons in the IGM. The re- 
quirements on the sky coverage, angular resolution and in- 
strumental sensitivity are demanding if a kSZ detection of 
the missing baryons is to be provided. Since the OV is mostly 
generated at moderately high redshifts, so if the cluster sam- 
ple does not extend deep enough, then the S/N drops accord- 
ingly (different survey depths are considered in the different 
panels of Figure (|4])). Note that we do not expect to find 
massive clusters at high (say z > 4) redshifts. 

If kSZ measure ment in halos are no t available, it has 
been suggested by IHo. Dedeo fc Spergell (120091 ') that kSZ 
templates may be built from the large scale matter distri- 
bution. Present reconstr uction algorithms like e.g., ARGO 
(IKitaura fc EnfilinI [iooj ) seem to reconstruct peculiar ve- 
locities accurately in scales of interest. This velocity field 
reconstruction could be an additional way to obtain the 
kSZ templates, and to detect the kSZ in clusters. In this 
case the S/N is much higher (dashed lines), and there is no 
need to have deep surveys: most of the S/N of the angu- 
lar cross-correlation is coming from the Poisson term. This 
term does not depend critically on depth, but simply relies 
on the discrete nature of halos. However, since the kSZ tem- 
plate would be a prediction of the model based upon the 
large scale structure observations, this exercise would still 
provide a test for the cosmological theory of bulk flows. 



6 DISCUSSION 

Our formalism allows us to study consistently the OV/kSZ 
temperature anisotropics introduced by both the gas col- 
lapsed in halos and the gas present in the IGM, regardless 
of the redshift range of relevance in each case. Our computa- 
tions differ from previous ones in mainly two aspects: (i) ours 
are full sky calculations, where no flat sky/Limber approxi- 
mation have been invoked and therefore predictions apply to 
all multipole range, and (ii) the kSZ computations take into 



account the fact that halos above a given mass threshold 
host only a fraction of the total mass, and hence the am- 
plitude of the kSZ they give rise to is limited. Indeed, even 
after accounting for the enhanced clustering of the massive 
halo population, we find that the kSZ generated by groups 
and clusters is below the OV anisotropies produced during 
and shortly after reionization. Only the Poisson term in- 
duced by the halo population reaches an amplitude that is 
comparable to the OV vv — dd term, although it should be 
possible to partially remove it by masking out the most mas- 
sive halos. By excising all halos above 2 x 10^'^h~^ Mq, the 
Poisson term should drop by a -^40 % (in the linear units 
of Figure ((3}). Distinguishing these two contributions may 
be a relevant issue, since the kSZ is mostly generated at low 
redshifts {z < 3 — 4), whereas the OV is one of the very few 
existing probes of the reionization epoch {z G [5, 12]). This 
gives significance to the cross-over of the total OV power 
(red filled circles in the left panel of Figure [3| above the in- 
trinsic linear CMB power spectrum at i ~ 4, 000. We remark 
that this OV computation is merely second order perturba- 
tion theory, and since the linear level of perturbations has 
been normalized via, e.g. WMAP observations, there are no 
free parameters for our (second order) OV predictions. Nev- 
ertheless, our computations do not observe patchy reioniza- 
tion, that is, the extra amount of anisotropy introdu ced by 
a non fully ionized IGM at 2: > 10. IZahn et all ([2005') claim 
that this phase of reionization should introduce fluctuations 
comparable in amplitude to the OV, and this would provide 
further relevance to the study of the anisotropies in these 
angular scales. 

This work is also motivated by the search of the missing 
baryons via their peculiar motions with respect to the CMB. 
If predictions from linear theory are correct, the baryons 
should comoving on scales of 20 - 40 /i~^Mpc, and hence 
the kSZ pattern generated in clusters and groups of galaxies 
should be very similar to that generated by the surround- 
ing gas, and can be used to probe the latter. However, it 
is not clear yet to what extent kSZ measurements in future 
CMB surveys like ACT or SPT will be contaminated by the 
intrinsic CMB and point source emission. In the ideal sce- 
nario in which errors in the kSZ estimates are unimportant 
for all clusters more massive than 2 x 10^*/i~^AfQ, then it is 
clear that the signature of missing baryons can be unveiled 
if the linear theory predictions on bulk flows are correct. 
One can consider to combine them with peculiar velocity 
held reconstructions obtained from LSS surveys, in which 
accuracy would be required at relatively large scales (10- 
20 /i~^Mpc). It should be noted that estimates of the kSZ 
would be needed for each bulk flow, and this can be achieved 
by combining kSZ/peculiar velocity estimates from different 
halos belonging to the same filament /supercluster. 

Our predictions have larger uncertainties in the kSZ 
case than in the OV case. We have assumed a Sheth-Tormen 
mass function for the halo abundance versus mass and red- 
shift. We also model their physical size and their gas profile, 
and more importantly, we have adopted the Mo & White 
prescription for their clustering bias. When comparing our 
results with previous w orks, we find t hat o ur kSZ power 
spectra are at the level of lSchafer et al.l (|2006l ) both in shape 
and amplitude: when imposing the same mass threshold 
(Af > 5 X 10^"^ h~^)MQ, our kSZ power spectrum is within a 
factor of 2, despite our different choice for erg, (0.817 versus 
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0.9). Further, we are introducing an effective bias for the 
halo velocities of 1.3, while it is not clear to us how close are 
halo peculiar velocities compared to li near theory expecta, - 
tions in the simulation they used. As ISchafer et al.l (|2006l ) 
pointed out, there may be a significant amount of extra kSZ 
power at smaller scales due to the pre sence of sub-structure 
within the halos, (|Springel et al.|[200ll ). This could partially 
be avoided by smoothing the signal within the halo solid 
angle, although part of this excess could be associated to a 
thermal (non-linear) component of halo velocities. The ad- 
dition of all this extra power could have an impact on the 
S/N estimations of Figure @, which increases significantly 
at high l-s. A lower mass for the kSZ estimation would then 
be required in order to achieve the same S/N limit. How- 
ever, despite all these sources of uncertainty, we conclude 
that future kSZ and LSS surveys would shed some light in 
the problematic of missing baryons and bulk flows in the 
late universe. 



7 CONCLUSIONS 

We have revisited the problem of the secondary anisotropies 
in the CMB introduced by the peculiar motion of baryons 
during and after reionization. We have addressed this prob- 
lem in the context of the search for the missing baryons 
and the measurements of bulk flows at low redshifts. We 
have paid particular attention in distinguishing the differ- 
ent sources of anisotropy, according to their redshift, in an 
attempt to isolate the kSZ/OV generated by gas in fllaments 
and superclusters at low redshift from its counterpart arising 
during reionization. For this purpose, we present the flrst all 
sky projection of the peculiar momentum of baryons from 
reionization to the present epoch. We consider two differ- 
ent scenarios: the OV generated by a smooth distribution of 
all electrons, and the kSZ generated by only those baryons 
located in collapsed structures. 

The former provides the largest amplitude, and con- 
tributes mostly during reionization {z > 6). Since it cor- 
responds to second order in perturbation theory, WMAP 
observations at the linear level already provide a normal- 
ization. Once those observations fix the set of cosmological 
parameters, the predictions at second order should also be 
fixed: the OV should cross over the intrinsic linear CMB 
temperature angular power spectrum at Z ~ 4, 000, with an 
amplitude close to 1.5 /iK (for a sudden reionization sce- 
nario at Zreio = 10.8 and r — 0.084). The amplitude of the 
OV should be above the angular anisotropy level induced 
by the kSZ in the halo population. The discrete character of 
halos (and the Poisson statistics associated to it) gives rise 
to angular fiuctuations that constitute the main contami- 
nant (in terms of kSZ/OV effects) to the signal generated 
during reionization. By masking out all halos more massive 
than 2 x 10^^ h~^MQ, this Poisson term should drop to a 
level of ~ 30% of the OV amplitude, leaving room for this 
window to the reionization epoch. 

The terms accounting for the velocity - velocity and ve- 
locity - density correlation of the halo population are well 
below the Poisson-induced anisotropy level, but are relevant 
since those terms express the correlation of velocities of halos 
with the velocities of the surrounding gas. We propose using 
kSZ estimates in halos (to be obtained by ACT or STP-type 



experiments) in order to build peculiar velocity templates to 
be cross-correlated to future CMB data. If the gas surround- 
ing halos is co-moving with them as the theory predicts, 
then a signal should arise from this test even after masking 
the contribution from the halos on the CMB maps. An ACT 
or SPT-like CMB experiment covering 10,000 (2,000) square 
degrees should achieve a S/N of 5 (2.5). The detection of this 
cross-correlation would provide evidence for the presence of 
missing baryons in the late universe and confirm the predic- 
tions of the theory of cosmological bulk flows. In case there 
are no quality kSZ estimates at the halo an gular positions, 
then, following IHo. Dedeo fc Spergell (|2009D . it would also 
be possible to build peculiar velocity templates from future 
deep LSS surveys. These would provide an estimate of the 
density held that can be inverted (within our cosmological 
frame) into a peculiar velocity template. By cross-correlating 
this template with future CMB maps one should be able to 
detect the kSZ and bulk flows at high signiflcance, if our 
understanding of cosmological bulk flows is correct. In sum- 
mary, cross-correlation studies of future CMB and kSZ/LSS 
data becoming a promising tool in the characterization of 
bulk flows and missing baryons in the low redshift universe. 
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APPENDIX A: THE 2-POINT ANGULAR OV/KSZ CORRELATION FUNCTION 

This Appendix outlines the computation of the fuU sky angular power spectrum of the projected peculiar momentum of 
baryons. It completes the appendixes B and C of HVJS in the sense that it considers the full peculiar velocity tensor 
(vjj i^q j) (and not only its diagonal part, as Equation Bl of HVJS implies). 

We follow the approach of HVJS for integrating kSZ temperature anisotropies along the line of sight. For a given population 
of halos, the kSZ introduced along the direction n reads 



^(n) = f dra(r)aT^W,a.(r-x,)ne,, (-^) = / dr a{r)aT f dx dM Vy,..(r-x)n.(M, r) UXi^j ^(x).(Al) 



The sum over j indicates sum over halos whose central electron density is given by riej and gas spatial profile by Wgasi^)- 
The coordinate r is in comoving units and can be regarded as our look back time coordinate and argument of the scale 
factor a(r) = 1/(1 + z{r)). The sum over j can be converted in an integral by introducing the halo mass function dn/dM. 
For notation simplicity, we shall implicitly integrate halo properties over the halo mass, dropping the mass dependence, and 
renaming the halo number density as nh(x). This transforms the integral in a convolution of the quantity 

C(x) = n,(x) = n,(x)/3(x) (A2) 



c 

with the profile Wgaa['x): 

^(n) ^ j dr a(r)aTn,,,(ri) {W^as^O W ^ j dr a(r)aTn,,,(ri) j W^^^^^^ e'^^ ^ (A3) 

The average central halo electron density is denoted by rie^h- The two point angular correlation function can then easily be 
written as 

C(m,n2) = j dr, a{r,)arn^Ari) J dr, a{r,)oTn^,u{r,) (^(^ ,-»k,.r, + .k..r. W^^^^^, W^;.,k. (^k^Ck,). (A4) 

The ensemble can be expressed in terms of the Fourier modes of the density and peculiar velocity field. Provided that ^(x) is 
the product of two quantities in real space, its Fourier counterpart equals 

^k = / (^/5q"h,k-q- (A5) 
Therefore, the ensemble average ((...)) in Equation (|A4|l becomes 

<?k,ek,) = j (§^3 (%«..k,-q/4.<.k.-q,)- (A6) 

One must address separately the different contributions on which this power spectrum can be decomposed, and this is the 
subject of the rest of the Appendix. 

In the rest of the Appendix, computations will be referred to a halo population. However, these same computations may 
be referred to a continuous and smooth distribution of electrons if certain changes are introduced in the equations. These 
changes involve: (i) the gas profile functions W ^^y^, the Fourier window functions associated to the peculiar velocity and 
density of halos {W^ j^i k' ^® introduced below) and the central halo electron densities n^^h must be dropped; (ii) the 
halo number density must be substituted by the electron number density field j^.; and (Hi) the halo clustering and 

velocity biases (6, 6„ , to be introduced below) must be set to unity. This continuous description of the electron field should 
be accurate when estimating the projection of the peculiar momentum of all baryons: the amount of collapsed baryons at 
high redshift is a very small fraction of the total, and even at present, less than a quarter of the total amount of baryons are 
found in collapsed halos more massive than 2 x 10^* h~^MQ. Regardless the amount of collapsed baryons, at scales larger 
than ~ 10 — 20 /i~^Mpc the matter density should closely follow the statistics of a moderately in-homogeneous Gaussian field. 



Al The vv term 

Let us first consider the case where the electron density or the halo number density are constant. This considerably simplifies 
Equation (|A6|) . since k-q ~ ?^h<5^(k — q) (with 5^ denoting the Dirac delta): 

(Ck.^k,)"" - = (2-)'5^(ki - k.)0.,,C„..:^(k, . nO(ki • n^) W^^^^^^ir,)W;^^^^^{r,). (A7) 

The factors denote the peculiar velocity time dependent growth factors along the the i-th line of sight. We are using the 
vorticity free linear theory expression for the Fourier modes of the peculiar velocity, vj^ — jD„k k^'k./^ with k the unit 
vector along k. We assign to velocities a Fourier window function k that smooths the contribution from scales smaller 
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than a halo if we are referring to a halo's peculiar velocity. In practice we set it practically equal to the gas window function 
W jj, which in general depends on the halo's mass and age. The particular choice of the velocity window function is not 
critical, since most of the anisotropy is coming from large scales. The matter power spectrum at present is denoted by Pm(k). 
After introducing this in Equation (IA4[) . one obtains 



C(ni,n2)= / dr-i a{ri)aTn^^h{ri)nh.i / dr-z a{r2)oTne,h{r2) nh.2 I P^.iPi;,2 e '^'''"^ ra) -PmW ^ 



(27r)3 fc2 

W^«.s,k('^i)W^;„.,k('^2) M/„,k('-i)Kk(^2)(Ai ■ k)(ri2 • k), (A8) 
which, after expanding the plane wave on spherical Bessel functions 

^-^k■r ^ ^(_i)'(2/ + i)P,(k . n)i,(fcr), (A9) 
and defining fi = a{ri)aTni reads like 

C(ni,n2)=^y" J driTinh,iV^,iW^^^]^{ri)W^]^{ri)ji{kri) J dr2i2nh.2V.u,2W*^^^^{r2)W*]^{r2)ji'{kr2)x 

^^^Pi(nr ■ k)(ni • k) Pi,{n2 ■ k)(n2 • k) {21 + l)(2l' + (AlO) 

We next define the line of sight integrals 

la = [ drd^flh.^V^.^W^,,^■^^{u)W^ ■^^{u)jl{kr,), (All) 



and apply the recurrence relation 

(a^) = m-i{f^) + + m+iM) ■ (A12) 

Finally, we perform the integral over the angles of k. We choose fii as the polar axis for k, so we need to remove the dependence 
of Legendre polynomials on n2 • k. We use the Legendre function addition theorem, that states that 



P,(n2 ■ k) = P,(k ■ ni)P(ni • 62) + 2 V ""j' cosm<?ij. Pr(k ■ ni)Pr(ni ■ 112). (A13) 

^ — ' (t + m)\ 

m — 1 

The azimuthal angle of k (here denoted as (j>]^) does not appear anywhere else in the integral, and therefore none of the terms 
in the sum proportional to cos m(f)y^ contribute to it. These manipulations yield 



I 



(^ 1) r-i-VV r-7-VV (^ 1)^ r-i-VV r-i-VV (^ ~t~ 1)^ ^VV ■-T-VV I ^ ■-j-VV r-7-VV 1 /' A 1 /I \ 

(27TTF " W+W " (27+ip • W+W ■ -^J ■ ^ ' 

Note that in the limit of high I and Ji '^ = I12 , the angular power spectrum multipoles vanish: 

lim Cr = Udkk^^ X - i - i + 0. (A15) 



i^oo ' TV J fc2 ' ' ' 4 4 4^ 

For this reason it will be required to compute the integrals T"^ with high accuracy, or the C^^-s will not drop to zero at high 
l-s as this limit requires. 



A2 The Poisson term 

This term observes the Poisson statistics of the number of halos along the line of sight. In this case, the power spectrum of 
reads like 

(^k.Ck,)" - / (§3 K,-q,4..q,)"(%/^q.) = (2-)' ^"(ki - - k2 + q2)n^Pf.(|k, - q,|) x 

^.,ki-q,W^;,k,-q, (2^)''5''(qi - qa) 0.,iI5.,2:^^(qi ■ ni)(q, ■ fia) W^^.q, VK;,q^ . (A16) 

The Poisson halo power spectrum is given by Phh(k) = 1/w, i.e., the inverse of the average halo number density. The Fourier 
density window functions ^ again limit the contribution of the halo power spectrum from scales much smaller than the 
halo size. They were approximated by Gaussian of characteristic size the halo virial radius. With this, the angular correlation 
function becomes 
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C(ni,n2) 



dri Ti / dr2 T2 nh-_ 



(2^) 



Let us next focus on the integral over q^, which shaU be denoted by K,: 



(A18) 



Let us take the polar axis for qj^ to be along k. In order to perform the angular integrals of the dot products present in 
Equation (|A18|) . we have to switch these products (q^ ■ fii), (q^ ■ 112) for others of the type (q^ • k) times other dot products 
independent of q. For that, we again make use of the Legendre function addition theorem: 

(qi-ni) = (qi-k)(ni-k)+cos0iP/(qi-k)Pii(ni-k), (A19) 
(qi • na) = (qi ■ k)(n2 ■ k) + cos02Pi'(qi ■ k)P/(n2 ■ k). (A20) 
The angles 01,2 refer to the angles between the planes containing the vector pairs (k,qi) and (k,ni), and (k,qi) and (k,n2), 

respectively. In this reference system where the polar axis is given by k, the vectors qj^, fii and n2 have each an azimuthal 
angle on the plane normal to k. We shall take the origin for this angle such that the azimuthal angle of fii is zero, and write 
02 = <?!>i + (5021 • According to the last two equations, the product (qj • ni)(qj • 62) reads 

2 r 



(qi • ni)(qi ■ n^) = (qi • k)'(ni • k)(n2 ■ k) + Pi'(qi ■ k) 



cos 0iPi (ill • k) cos 02^/(112 • k) 



cos0iP/(qi ■ k)Pi (fii • k)(qi • k)(fi2 • k) +cos02Pi (qi ■ k)Pi (fi2 • k)(qi • k)(fii • k) 



+ 



(A21) 



The last two terms proportional to cos 01,2 will give no contribution once the integration on the azimuthal angle of q (0i) 
is done, and therefore can be neglected. Since 02 = 0i + 502i, then cos 02 = cos0i cos502i — sin0i sin502i. The product of 
cos 01 cos 02 will therefore give rise to a term proportional to cos 0i sin 0i that again gives zero contribution to the integral 
over the azimuthal angle of (0i). Thus, effectively, we are left with the term proportional to cos'^ 0i: 



(Qi ■ ni)(qi ■ "2) = (cii • k)^(fii • k)(fi2 •!<;)+ P/(qi ■ k)cos0i 



cos 5021 Pi (fii ■ k)Pi (fi2 • k) 



(A22) 



We again invoke the Legendre function addition theorem on the triad formed by the vectors k, fii and fi2 in order to manipulate 
the expression within the square brackets of the last term on the rhs of the last equation. It now becomes 



(Qi • fii)(qi ■ 112) = (cii • k) (fii • k)(fi2 • k) + Pi (qi • k) cos0i 



(fii ■ 112) - (fii ■ k)(fi2 • k) 



(A23) 



This expression can be plugged into Equation (|A18|I : 



qidqi 
iM' Jo 



^dM,i|W^.,qinM/,k-qj'^^ 



COS^ 0llPl'(^i,JP(fii ■ fi2) + 



[nl^ - cos^0ijPi (^iqJI^) (fii ■ k)(fi2 ■ k) 



(A24) 



where = qj^ ■ k. We next assume that, on average, halos are spherically symmetric, so that Wg = Ws,k and = W^^k- 

This simplifies the last equation. 



K. 



Jiik){h^ ■ fia) - J2(fc)(fii ■ k)(fi2 • k) 



with J'i{k) defined as 



and 1/2 (fc) as 



Pm(qi 



J2{k) 



(27r)3 



rfM,i \W..qA'\Wg,^_^/^'^ 



f4l - 1-Pl (Mqi)l^COS^01 



This also simplifies the expression for the angular correlation function (Equation IA17|) : 



C(fii, fi2) 



dri fi / dr2 T2 nh,i 2?„,iD„_2 



dk 



ik.(ri-r2) 



W 



Ji{k)(ni ■ n2) + (fii ■ k)(fi2 • k)j2(fc) 



(A25) 



(A26) 



(A27) 



.(A28) 
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Let us split it into two contributions, according to their angular dependence. The first contribution Cj(ni,n2) contains the 
term (fii • 62), the second one {Cglni, 0.2)) contains (fii ■ k)(n2 ■ k). The procedure now is very similar to that used for the 
vv term. We first expand the plane wave as in Equation (|A9|l , and define the line of sight integrals Xf^^ as 



T/;, = J dn h ^nh,r W^^^yP^,,ji{kr,). (A29) 

A2.1 The CA(ni,n2) part 
The contribution from Cji(ni,n2) then reads 

Cj(ni,n2) = + 1){21' + / -0- Mk) I^,!^^ Pi{Ai ■ k)P,,(n2 ■ k)(fti ■ 112). (A30) 

The addition theorem applied on Pi'(n2 • k) and the integral on the azimuthal angle (that removes terms of the type cosm^j^.) 
leave the last expression as 

Cl(ni,n2) = ^(2/ + 1)(2Z' + [ ^li^fa Pi{ni ■ k)P,'(fti • k)Pi,(ni • n2)(ni • 112). (A31) 

1,1' 

Finally, the recurrence relation (|A12|I yields the final expression for C^(ni,n2), after integrating in /ij^: 

C7l(ni,n2) = ^^P(fti-ft2)CrA = 5]^Pi(fti-ft2) (J) J fc'dfc ^i(fc) f ;jA4x.^i,iir+i,2 + TTr7T^r-i,iXr_i,2 Y(A32) 



In this case, this contribution does not vanish at high I 



21 + 1 2; + i 



lim C^A = - I k^dk Ji{k)\irf. (ASS) 



A2.2 The C'£ {ni,n2) part 
This contribution is written like 

Ci^(ni,n2) = -^(2« + 1)(2/' + l)(-i)'-'' / J2{k) ifilfa Pi{Ai ■ k)P'(n2 • k)(ni • k)(n2 • k), (AS4) 

1,1' •' ^ ' 

which, after using the Legendre polynomial recurrence relation (Equation IA12I) and integrating on the angles, yields 
C£(ni,n2) = -^^^PKni-n2)Q^B = -^^^P(ni-n2) (J) f k^dk J2{k)x 

As it was the case for the vv term (which showed an identical dependence on fii, 62), this contribution will vanish at high I 



^P 2 J, , 2 ^ / , \ /'tP\2 /I 1 1.1 



fim C[^B = - dkk' J2{k) (irr X t-T-T + T ^O- (A36) 
i—oo TV J \4444' 

Therefore, we find that limz^oo Cf = Cf j^. 



We have concluded the computation of the linear terms, and we next address the computation of the second-order terms. 
In this Section, the electron density will be assumed to trace the dark matter density field. This means that, when looking 
at the kSZ generated in halos, the number density of halos will trace the local, large scale, dark matter density field, in such 
a way that fluctuations in the halo number density will be a bi ased mapping of the fluctuation of dark matter density (and 
this bias factor was approximated by that of IMo fc White! (|l996t )). I.e., Jt-^^ k = bT)sS■^^W^ j^, with b the bias factor, the k 
Fourier mode for the halo number density distribution, Jj^ the Fourier mode of the dark matter density contrast at present, 
and T>s its linear growth factor. When considering the OV generated by the IGM on the large scales, this assumption is 
expressed as = fie{r)T>s 5-^ with ne(r) the background electron number density. 
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A3 The vidi — 1)2^2 term 

This term observes the contribution of the configuration 



(A37) 



As we show next, this contribution vanishes. For a real field A(x), we have that = ^* j^, so the last equation can be 
rewritten as 



(2.r)«5^(kx)^P™(gi)(q, ■ m) 27r)^5°(k2)^P„(g2)(q2 • A,). 



If we choose fii as the polar angle for q^ {i = 1, 2), then Hi = (q^ • hi) and we end up with an integral of the type 



qfdq^ Pm{q 



rf/ii ^J.i = 0. 



(A38) 



(A39) 



Therefore, this term provides no contribution and will be ignored hereafter. 



A4 The vv — dd term 

We next compute the contribution of the configuration 

= J (2^(2^ <^*li^q2)K.ki-q,^h.k,-q,) 



within the two-point angular correlation function 



(A40) 



The ensemble averages of Equation (|A40|I yield 



dc\-^ dqj 



(27r)3 (27r)3 
Pm(gi), 



(27r) 5 (q^ - q2)l'i.,il'i,,2 2 — M^^',ql^^h,lW^u,q2"h•2 (^li ' ni)(q2 ' "2) X 



(27r)3 (2^)3 

(2^)^5^(ki - qi - k2 + q2)%,i&i%,2&2P™(ki - qi)W^,,ki-q, W^;,k,-q, 
The presence of these Dirac deltas significantly simplifies Equation (IA41[) : 

rfQi Pm(gi), 



(A42) 



/ ^^|W^.,qj'(qi-fti)(qi-ft2)P™(k-q,)|W^,k-qj' 



(A43) 



There are functions inside the integral over qj^ which depend upon k — qj^ . Just as for the Poisson term, this suggests using k 
as the polar axis for the angular integration of q^, /i^j = k-qj^. Thus we have to apply the Legendre function addition theorem 
and follow the same procedure as in Equations (IA19I - |A23|) . After introducing the usual expansion of the plane waves, this 
gives 

C—''{h„h2)=J2^-^y-''J j^I--'Xr'"'Pl{i^i-i^)Pl'ii^2-k) J :5:i^|M/..q/P„(k-q,)lW^,_k-qj' >< 



qi • ki)^(ni • k)(n2 ■ k) + cos<j!>iPi (k • qi) I ni ■ n2 - (rii ■ k)(n2 • k) 



(2Z + 1)(2«' + 1), 



(A44) 



where the T,"""'*'' = J dri ji{kri)Tinh,iW^^^ i^biVy^iVs^i are the usual line of sight integrals. The integration on the azimuthal 
angle of q^ can be done trivially, but the integration on its polar angle must be done numerically. These angular integrations 
introduce two fc-dependent functions, A{k) and B{k): 

^vv — ddf- " \ \ ^ / — f ^k ^vv — dd'-rvv — ddrt/-x 1 \ rj i\ f Ql^qidflq^ Prn{qi) ^j-r^ 1 2 j-* /i „ Wxxr \~ 

C (ni,n2) = y (^^'.1 ^i'.2 -P!(ni-k)Pi,(n2-k) j—^^ _| V[/„^qJ P„(k - qi)|M/^_k_qJ 



(2^) 



(ni . k)(n2 • ic) ( (k ■ q.f - i (Pi^(k . qj)' ) + \ {pI{\^ ■ qj)' (ni • 112) 



(2/ + l)(2/' + l) 



dk 



^'^,(n,.k)P,(n2-k) y gr^li:^ilV„,qJ^P^(k-q,)|W',k-qJ 



qxdqx Pm{qi) 
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A{k){ni ■ k)(fi2 • k) + S(fc)(ni • na 



(2/ + l)(2Z' + l). 



(A45) 



These two functions will give rise to two different contributions for the vv — dd term, 112) and ''''(iii, n2)- 



A4.I The CI" -'^'^ {ill, (12) contribution 



As already done for the vv term and one of the two contributions for the Poisson term, we first apply the Legendre polynomial 
recurrence relation (Equation IA12[l in order to express all the dependence on k • ni and k • n2 as arguments of Legendre 
polynomials: 



/P,_i(ni-k) + (/+l)P,+i(nrk) 



/'P,,_i(n2-k) + 0' + l)Pr+i(ft2-k) 



A{k).{A4:Q) 



The next step uses the Legendre function addition theorem (Eauation lA13|) to rewrite Legendre polynomials whose argument 
is k • n2 in terms of Legendre functions of arguments k • fii and fii • n2. The integral on the azimuthal angle of k removes 
all terms proportional to cosm^j^, so in practice we simply substitute the Legendre polynomials Pi(k -112) by the product 
Pi(k • n2)P;(ni -112). The polar angular integral of k does the rest: 



Ca (ni,n2) = 2^ Pi (ni • n2)C,^ " ' "2 



- / k^dk A{k) 

TT 



2 



I 7orT~n2 '+i>2 ^ ^o; , i\2-'-!-i,2 I +-''1-1,1 I ^ /r,7 , iN^-'-i-i.a + ■ 



^vv — dd 
)'-^i-1.2 



.(A47) 



As for every term with the angular dependence of the type (fii • k)(n2 ■ k), its contribution vanishes in the limit of I ^ <x and 

■vv- 
1,2 



'-T-w — dd ^vv — dd 

-''1,1 = 



hm C^r"' = - I dkk" A{k) (i^^-^y xfi-i-i + i|=0. 



(A48) 



A4.2 r/ieC^"" (ni,n2) contribution 

In this case we again have to use the addition theorem of Equation (IA13|) to obtain 
Cl"-^''{ni,n2) = + + 1) J j0^^r''X'7'"' l3{k)^PAni ■ n2)(ni • n2), 



(A49) 



where represents the Kronecker delta = 1 if i = j and — otherwise). Finally, an application of the relation 



(|A12|) gives 

Cr-''^(ni,n2)=^^P(ni.n2)C,V'' = E^^'(^^-»2)^ / k'dkB{k) 
I I '' 

Hence, this term gives most of the contribution of (7™"'*'* at high l-s: 



y-rvv — dd'j-vv — dd ^ ~t~ 1 ^ ^vv — dd^vv — dd ^ 
-'■l + l,! -'■1 + 1,2 21 -\- 1 '-1.1 i-1,2 



lim C. 



= lim C," 



2; + 1 



(A51) 



.(A50) 



A5 The vd — vd term 

In this final section, we compute the contribution coming from the configuration 

(^kA7''^'^l (§y3(^ (/5qi";,k.-q,)K,k,-q/4.) 

within the two-point angular correlation function 



(A52) 



— — «-''^^"^+''^^"^W^^„.,k,W^;.,k.<Ck,Ck,)"^-''''.(A53) 



C (ni,n2)= / dn a[ri)aTn^^h[ri) J dr2 a(r2)aTne,h(r2) j^;^ j^;^ ' 

This time, the computation of the ensemble averages of Equation (|A52|) is slightly more complex. We first obtain, after 
canceling the integral on k2 by using one of the Dirac deltas. 



i^kAf"^' = / (^(2^)'5''(qi +ki -qi -k2)n,,,262l)M2'.,iW^.,qiW^;,k,-ki+q,(qi ' "i) 



Pm(gi) 
91 
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ki - qi 



Pm.(ki - 



(A54) 



Again we have an integral on qj of functions whose argument is ki — qj^ . This means that we have to translate all arguments of 
the type (q^ • ni) and (q^ • n.) onto the quantities (qj^ • ki) and (fii ■ n.) and choose ki as the polar axis for the q^ integration. 
By using Equations (|A19I |A20)) and neglecting the terms proportional to cos0i, cos(;62 (which give no contribution after the 
integration on the azimuthal angle of q^), one readily finds 

-Pm(gi) -Pm (kl - qJ 



91 



|ki - qil 



fci(ki •ni)(ki •n2)(qi -kO-gi {{q^ ■ ki)'(ki ■ ni)(ki ■ fia) +cos0i cos<^2|Pi'(qi ' kil'Pi (ki ■ ni)Pi'(ki • m)) 



.(A55) 



As we did for the Poisson term, we now express cos (j}2 in terms of and 5(j]2i = 02 — 01 and make use of the relation 
cos 5021^1 (kl ■ ni)Pi^(ki ■ £12) = (fii ■ 62) — (ki ■ ni)(k2 ■ n2)- After integrating in the azimuthal angle of q^, this leads to 



<?k,ek,>" = (2vr)'5°(ki-k2) 



qldqidfiq, a a -n -n -n -n w w* w w* Pm(gi) -Pm(ki - qJ 

■nh,inh,20ib2Vs^2Vv,iT>i^iVv,2Wv,ci^ Ws^d^ '^a.ki-qj ki-q^ ~ '• — ^ 



(27r) 



(2^) 



(kl ■ ni)(ki ■ 112) (^fciM,i -giMqi + y1P/(M9i' 



|lPi^(M,i)l'(ni-n2) 



gi Iki-qJ 
(A56) 



We plug this expression into Equation HA53P after expanding the plane wave onto spherical Bessel functions and defining 
Ti = a{ri)cjTne,h{ri) and the line of sight integral Xj"^"""* = ^ij''^'' = / dn ji{kri)TiW^^^ -^fih.ihiDv,i'Ds^i. We obtain 



/^vd — vd / " " \ ^ / — I n-vd—vdn-vd — vd 7-1/" \ n / T 

C (m.na) = 2^(-i) 72;p^^i.i- ^i',2 Pi(ni • k)P,/ (n2 • k 

1,1' •' 

(k- ni)(k- n2) {kl^lq^ - qifi]^ + ylPi (m9i)I^ 

EH)'"'' / 7^l":r"-P^(fii-k)P.(n2-k) (2^) 
, ,, J \ ) 



qldqidfiq^ P'm(gi) P,7i(k - qJ 



ylPi (M<7i)l^(ni ■ na) 



(27r)3 gi |k-q,|2 

(2; + l)(2;' + l) = 



(27r)x 



(k-ni)(k-n2)C(fc)-0(fc)(ni-n2) 



(2; + l)(2Z' + l). (A57) 



The two functions C{k) and ©(fc) introduce two different angular dependences of the correlation function ((k • ni)(k ■ 112) 
and (fii • 62), respectively), and give rise to two different contributions to the total vd — vd correlation function, that will be 
denoted by C^'^~'"*(ni, 62) and C]^'*~'"'(ni, 62), respectively. 



A5.1 The CI''""' {ill, n2) contribution 



As usual, we first use Equation (|A12p to remove the factors of the type xP; (a::) and then apply the Legendre function addition 
theorem to switch arguments of (k • 112) into arguments of (k • fii) and (fii -£12). All terms proportional to cosm^j^ give no 
contribution, so under the integral sign we are allowed to replace Pi(k ■ fi2) by Pi(k • fii)P;(fii • £12). This results in 



k^dkC{k) 



.-rvd—vd I (^ ^-pvd—vd ^(^ ^^ 1) ^vd—vd \ , ■-j-vd — vd 

-''i + l.l I /nl , i\2 i + 1,2 ^ 77rrTTv2-''t~l,2 ) +-^-(-1,1 



(2/ + 1)2 



(2Z + 1)2- 



And again, in the limit of high I and 



we have that 



+ 1) -rvd—vd . 
jM-\,2 + 



I' 



{2i + iy 



{21 + 1) 



7-^1-1,2 



.(A58) 



lim cr 



dk k" C{k) (j^"'*~"'')2 X 



(A59) 



A5.2 The Cl''"''^{ni,h2) contribution 

This case is again identical to the C'^~''''{ni,n2) contribution. The final result is 



Cr'"^(fii, n2) = J2 ^^'("1 • "2)C, 



■vd — vd 
l,B 



y ^^^'(ni • - / k'dk {-v{k)) 

l 



^vd — vd^vd — vd ^ + 1 . '-rvd — vd^vd — vd 
-^1 + 1,1 -^1 + 1,2 o7 I 1 ~'"-'';-l,l -^1-1,2 



I 



21 + 1 



21 + 1 



(A60) 
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This term gives most of the contribution of Cj"'' at high l-s, and it is negative. It is always below the vv — dd term, but it 
reduces the (positive) amplitude of that term. 

lim Cif" = - f dkk" i-Vik)) ix^^-^-'i'if = lim cf""^. (A61) 



